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Abstract 

Interaction  of  bianisotropic  particles  in  plane  regular  arrays  is  under  investigation  in  this 
paper.  We  give  a recipe  on  how  to  use  the  particle  polarizabilities  and  the  interaction 
constants  obtained  from  approximate  models  so  that  the  energy  balance  is  satisfied  and 
the  physics  of  the  phenomena  is  kept.  Numerical  examples  are  given  for  arrays  of  omega 
particles. 

1.  Introduction 

To  solve  a diffraction  problem  for  an  array  of  scattering  particles  one  should,  at  first,  know 
the  properties  of  an  isolated  inclusion  given  by  its  polarizability  and,  second,  how  the  inclu- 
sions interact  in  the  array.  The  polarizabilities  as  well  as  the  interaction  fields  often  cannot 
be  calculated  exactly.  The  aim  of  this  paper  is  to  give  a method  which  will  allow  us  to  use 
the  polarizabilities  and  the  interaction  coefficients  obtained  by  approximate  models  so  that  the 
reflection  and  transmission  coefficients  will  satisfy  the  energy  conservation  principle. 

2.  Energy  Conservation  in  Bianisotropic  Arrays 

Consider  a plane  regular  array  of  scattering  particles.  We  will  assume  that  the  particles  may  be 
represented  as  combinations  of  electric  and  magnetic  dipoles.  Thus,  every  particle  is  character- 
ized by  its  dyadic  polarizability  factors  E: 

P = aee  ’ Eloc  + aem  " H-ioc  . . 

m — ame  ' Eioc  + OLmm  ' H(oc 

We  assume  that  the  array  is  excited  by  a normally  incident  plane  wave  with  the  fields  Eext  and 
Hgxt-  Every  particle  is  excited  by  the  local  fields 

Eloc  = Eext  + Pe  • P>  Hloc  = Hgxt  + Pm  • HI  (2) 

Here  (3e  and  /3m  are  the  interaction  dyadics.  These  dyadics  take  into  account  interaction  of  the 
particles  in  the  array.  Under  our  assumptions  the  array  of  electric  dipoles  does  not  produce 
any  magnetic  interaction  field  and  the  array  of  magnetic  dipoles  does  not  produce  any  electric 
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interaction  field.  Due  to  this  there  are  no  cross  terms  in  (2).  As  it  was  shown  in  our  recent  work 
[1]  the  interaction  dyadics  for  the  considered  problem  can  be  represented  as 


fie  = Re(/3e)  + j 


W qHquc 
6tt 


■quj  = 

3 2 SQIu 


fit 


m 


U)  as 

It 


2rjSo 


(3) 


Here  Sq  is  unit  cell  area.  As  one  can  see,  the  last  terms  in  the  above  relations  correspond  to 
the  plane  wave  field  contribution.  We  can  express  the  local  fields  in  terms  of  the  induced  dipole 
moments: 
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The  external  fields,  as  follows  from  (2),  can  be  written  as 


(4) 


Eext 

Hext 


— (&ee  ~ 


= ( 


amm  ame  * aee  ' * (m  ame  * &ee  ' P)  fim  ’ m 


®em  ' amm  ' ® me)  * (P  &em  ' ®mm  ‘ m)  fie  ' P 


(5) 


The  total  averaged  fields  (plane  wave  fields)  in  the  array  plane  read 

Etot  = Eext 


Htot  — Hext  j2r)s0Vm-  Z 


' P — Zee  ■ J + Zem  • J rn 
u3-r  - m 7me  ‘ J "h  Zrnrn  • Jm 


(6) 


Here,  the  currents  J and  Jm  axe  not  the  average  surface  electric  and  magnetic  currents.  These 
vectors  can  be  arbitrarily  directed  and  they  represent  the  normalized  electric  and  magnetic 
dipole  moments.  The  dyadic  coefficients  in  (6)  can  be  easily  identified  from  the  above  formulas: 
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Although  different  terms  have  different  dimensions,  we  use  the  same  notation  Z for  all  of  them. 
Indeed,  only  Zee  has  the  meaning  of  impedance. 

Let  us  now  suppose  that  the  particles  have  no  dissipation  losses.  Then,  the  energy  conser- 
vation condition 

Re{Etot  ■ J 4-  Jm  ■ Htot}  = 0 


can  be  written  in  the  dyadic  form  as 


Zee  + Ze 
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Since  this  should  be  valid  identically  for  all  currents,  we  find  that  in  lossless  arrays 


Zee  + Zte  — 0) 
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where  f denotes  the  Hermite  conjugate.  Consider  reciprocal  particles.  Then  Zee  and  Zmm 
are  symmetric  dyadics.  Thus,  (13)  means  that  these  dyadics  are  purely  imaginary  (and  the 
bracketed  expressions  are  purely  real).  In  other  words, 
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and  similarly  from  (14): 

(17) 

(18) 

The  last  two  relations  are  equivalent  since  in  reciprocal  media  Zme  — —Zern.  For  a special  case 
of  omega  particles  with 

See  = flee  X0X0  + a^YoYO,  Smm  = a,mmZ0Z0,  Sme  = ameZ  0y0,  Sem  = “OmeY  0*0  (19) 

relations  (15)  and  (16)  give 
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Conditions  (17)  and  (18)  lead  to  the  same  result,  which  reads 

Re{  yy  2 }=0 

l + afne  J 


(21) 


3.  Reflection  and  Transmission  Coefficients 

Let  us  now  make  the  final  step.  Our  mean  goal  is  to  find  the  reflection  and  transmission 
coefficients.  It  is  easy  to  see  that  they  can  be  found  in  terms  of  the  introduced  parameters 
(7)-(10).  The  parameters  are  expressed  via  the  particle  polarizabilities.  In  the  usual  practice 
we  have  the  polarizabilities  dyadics  found  from  the  antenna  model  or  numerically.  This  gives 
approximate  results  and  the  values  of  alphas  are  not  quite  correct.  We  should  somehow  “correct” 
the  polarizabilities  to  find  the  reflection  and  transmission  coefficients  which  satisfy  the  energy 
conservation  law. 

From  the  other  hand,  the  correction  should  not  lead  to  significant  difference  in  the  array 
reflective  properties.  If  we  simply  skip  the  scattering  terms  from  the  polarizabilities  and  the 
interaction  dyadics,  i.e.  we  take  aee  and  amm  as  purely  real,  aem  and  ame  as  purely  imaginary, 
and  leave  only  plane  wave  contribution  in  the  imaginary  parts  of  betas,  we  easily  satisfy  the 
energy  conservation  conditions  (15)— (18).  But  the  frequency  behaviour  of  the  array  reflection 
will  change  dramatically.  It  follows  from  (7)-(8):  equating  the  scattering  terms  in  the  alphas 
to  zero  we  change  the  values  of  dyadics  Z.  It  suggests  us  to  apply  additional  conditions  to 
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avoid  that.  These  conditions  will  keep  the  non-scattering  terms  of  (7)-(8)  unchanged  when  the 
scattering  terms  in  alphas  and  betas  are  skipped.  Doing  so  we  obtain  a system  of  equations  on 
the  new  “corrected”  polarizabilities  aee,  amm,  aem,  and  ame.  In  more  detail,  the  system  is 


M«ee)  = Im(5^nm)  = 0,  Refa^)  = 0 


(22) 
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The  “corrected”  polarizabilities  together  with  the  “corrected”  betas  (only  plane  wave  contri- 
bution is  included)  will  give  the  reflection  and  transmission  satisfying  the  energy  conservation 
law. 


4.  Numerical  Results 


We  have  numerically  investigated  the  case  of  a double  array  of  omega  particles  under  the  plane 
wave  excitation.  The  polarizabilities  were  obtained  from  the  antenna  model  described  in  [2]. 
The  reflection  coefficient  values  via  the  frequency  together  with  the  energy  balance  plot  are 
given  on  Figure  1 and  Figure  2.  The  solid  lines  correspond  to  the  corrected  alphas  and  betas, 
the  dashed  lines  correspond  to  the  original  ones. 


Figure  1:  Reflection  coefficient  as  function  of  Figure  2:  Energy  conservation  as  function  of 
frequency.  frequency. 
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